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مقدمه 

فا کون ووی‌هان بای اب سیون سا مسا یت 
سازی پيشنهاد شده و بکار گرفته شده است. برخی از اين 
روشهای حل مسائل بهینه سازی توانائی رسیدن به 
جوابهای بهینه مطلق را دارا می باشند و برخحی دیگر مانند 
الگو ریتم های فراکاوشی (عصطاتتمعله عتاعتسعط ۸۵۲2) 
تنها به رسیدن به جواب‌های خوب و مناسب بسنده 
کف فایتفه فدانش وهای که انا دنق با 
جوابهای نزدیک به جواب بهینه را دارا می باشند را می 
تفت اه وافست خاسست: که برشی از سای بر 
سازی شدیدا غیرحطی می باشند به گونه ای که 
الگوریتمی صریح برای حل آنها وجود ندارد [1]. از 
طرفی در بسیاری از مسائل واقعی. از جمله مسائل 
مهندسی آب تصمیم گيرنده ها تنها نیازمند جوابهای 
خوب و مناسب (نه لزوما جواب بهینه مطلی) می باشند. 
بر این اساس امروزه استفاده از روشهای کاوشی و یا 
الگو رِ یتم‌های تکاملی (قصطاتمعاه ممماانان ۲ که 
تضمین کننده رسیدن به جوابهای نزدیک به جواب بهینه 
می باشند. مورد توجه قرار گرفته است . 

الگوریتم های فراکاوشی نامی معمول برای 
روش‌هایی است که عموماً برگرفته از طبیعت هستند و 
اساسا سعی در ترکیب اصول اولیه روشهای کاوشی برای 
رسیدن به یک جستجوی موثر و کارا در محدوده مورد 
نظر دارند. در سالهای اخیر الگوریتم‌های فراکاوشی با 
پیشرفتی سریع در حل مسائل پیچیده و مشکل بهینه 
سازی مورد استفاده قرار گرفته اند. الگوریتم نورد شبیه 
سازی شده (8۵) (عصتحمعصمم عطاق الگو ریتم 
ژنتیک (0۸۵) (صطاندهع۸1 ع1اعظ06). جستجوی ممنوعه 
(15) (62700 1800 جستجوی موضعی اکتشافی 
(صاظ) رطعم لقهما ۵تامتمام۲) و الگو ریتم بهینه 
سازی جامعه مورجه ها (۸00) ( بعمامن اصظ 
صطاتتمعله حمتاعمتا00) را می توان به عنوان نمونه 
هایی از اين الگوریتم ها نام برد. الگوریتم بهینه سازی 
جامعه مورجه ها اولین بار توسط 001001 و همکارانش 
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در سال ۱۹۹۱ پيشنهاد شد [2]. الگوریتم سیستم مورچه‌ها 
(5ه) (صعادر 5‏ ۰ ۸۶) اول‌ین الگ‌وریتم از دسته 
الگوریتم‌های پیشنهادی جامعه مورچه ها بود. در طول 
سالیان اخیر الگوریتم های مختلف دیگری نیز بر پایه اين 
الگوریتم پایه اولیه ارائه شده است که در ادامه به برخی 
از آنها اشاره می شود. از جمله» ۲0۲۱۵0 و 0202706112 
در سال ۱۹۹۷ الگوریتم سیستم جامعه مورچه ها (۸09) 
(عاولرو. . بزجمآم ۰ ۸۵۲) را معرفی نمودند [3]. 
عصزه130110 و همک‌ارانش در سال ۵۹ با اصلاح 
الگوریتم سیستم مورچه هاء الگوریتم سیستم مورچه های 
ترتیبی (۸9۵1) (صعادر٩‏ ۸ ۵۵60 ) را پيشنهاد 
نمودند [4]. 5001216 و ۲1005 در سال ۰۲۰۰۰ الگوریتم 
سیستم مورچه هىای بیشینه - کمینه (۸۹) (-۱۸۵ 
«عاو5 اه۸ ۷2) را معرفی نمودند [5]. از جمله دیگر 
الگوریتم های پیشنهادی. می توان به الگوریتم سیستم 
مورجه های بهترین -بدترین ( ۸0 ۳65۷۷0۲5۹۲ 
8 [6]؛ الگوریتم سیستم مورچه ترکیسی ۲۱۷9714 
مدرگ نصص [7] و الگوریتم مقید جامعه مورحه ها 
(صطاتمعاه ممتاععنصتاع0 وجمامت امه 0عصتحعومن) 
[8] اشاره نمود. در ادامه به برخحی از کاربرد های الگوریتم 
بهینه سازی جامعه مورچه ها در حوزه مهندسی آب از 
زمان پیدایش این الگوریتم اشاره شده است. اولین کار 
گزارش شده دراین حوزه مربوط به ۸۰0۵5۲۷0۲ و 
هار انش هر سا ۱۲۰۷۱ [9] مس نت که از انش 
الگوریتم برای تخمین پارامترهای هیدرولیکی خاکهای 
غیراشباع استفاده نموده اند. پس از آن 51000500 و 
همکتاراین فرستال ۱۳۸۰۲ 101 با انتفاهه از ارسته 
الگوريتی مساله طراحی بهینه شبکه آب را مدل نموده و 
پارامتر های آن را محاسبه نمودند.۱]2167 و همکارانش 
در سال ۲۰۰۳ [11] ازاین الگوریتم به منظور طراحی بهینه 
یک سیستم شبکه توزیع آب استفاده کردند. آنها با مقایسه 
نتایج به اين نتیجه رسیدند که این الگوریتم یک الگوریتم 
قانل رقایته با الگرونش وکا در شاه ظراخی نیت 
شبکه توزیع آب می باشد. 760010 و همکارانش در 


سال بیست و سوم شماره یک. ۱۳۹۰ 


محمدهادی افشار- ابراهیم رضایی سنگدهی- رامتین معینی 


سالهای ۲۰۰۳ و ۲۰۰۵ [12,13] ابتضدا از الگسوریتم 
مورچه‌های بیشینه - کمینه در مساله طراحی بهینه شبکه 
آب استفاده نمودند و سپس عملکرد الگوریتم سیستم 
مورچه ها را در مقایسه با الگوریتم مورچه های بیشینه - 
کمینه در مساله طراحی بهینه شبکه آب مقایسه نمودند. 
افشار در سال ۲۰۰5 ابتدا از الگوریتم مورچه های 
بيشینه - کمینه در مساله طراحی بهینه شبکه آب امستفاده 
نمود [1] و سپس با بکار گیری مکانیزم تظریف تطبیقی 
قطعی (۵۲) (رتصعصم‌صی؟ ع۶تامعک متاوتمتماعر) 
در الگوریتم مورچه هاء عملکرد این الگوریتم را بهبود 
بخشید و اثر این مکانیزم را در طراحی بهینه سیستم 
کنترل سیلاب مورد بررسی قرار داد[14] . همچنین افشار 
و معینی در سال ۲۰۰۸ [8] الگوریتم مقید جامعه 
مورچگان را پيشنهاد نمودند و کاربرد آنرا در بهره برداری 
بهینه از مخازن سدها مورد بررسی فرار دادند. 

الگوریتم استاندارد اولیه پیشنهادی جامعه مورچه ها 
اساسا برای حل مسائل بهینه سازی گسسته, که در آنها 
متغیر های تصمیم تنها مقادیرخحاصی را می پذیرند» ارائه 
شله: ابیت بتایراین به .کار گیری این الگوریتم پبرای حل 
مسائل پیوسته نیازمند تبدیل فضای پیوسته جستجوی 
مساله به فضای جستجوی گسسته از طریق گسسته سازی 
متغیرهای تصمیم پیومته مساله میباشد: کقیست رات 
نهایی و عملکرد مکانیزم جستجو تا اندازه زیادی به ابعاد 
بازه های گسسته سازی متغیرهای تصمیم پیوسته مساله 
وانسته: اسیت؛ ور کته سازی درشت متغیرهای تصمیم 
مساله ممکن است گسسته سازی ها به گونه ای انجام 
شده باشند که شامل نقاط بهینه مساله نباشند و گسسته 
سازی های ریز موجب افزایش تلاش محاسباتی و بعضا 
کاهش کیفیت جواب می گردند. مکانیزم تظریف تطبیقی 
قطعی پيشنهاد شده در این مقال» جایگزین بسیارخویی 
برای جلوگیری از گسسته سازی های زیاد متفیرهای 
تصمیم مساله به منظور ساختن فضای جستجو مساله می 
باشد. زیرا این روش نیاز به هزینه عملیاتی کمتری دارد و 
به جوابهای دقیق تری نیز دست می‌یابد . 
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در اين مقاله کاربرد مک‌انیزم پیشنهادی (تظریف 
تطبیقی قطعی) در الگوریتم سیستم مورچه های بیشینه - 
کمینه به منظور بهینه سازی توابع ریاضی و مسائل بهره 
برداری بهینه از مخزن مورد بررسی قرار گرفته است. 
نتایج بدست آمده بانتایج پدستت امنلده ای بکناو گیری 
روشهای دیگر و با نتایج بهینه بدست آمده از مدل سازی 
مساله در نرم افزار 11300 (نسخه ۸) مقایسه شده است. 
مقایسه نتایج نشان دهنده آنست که مکانیزم پیشنهادی از 
قابلیت بالایی در دستیابی به جواب بهینه یا نزدیک به 
بهینه با تلاش محاسباتی کم برخوردار بوده و می تواند 
در حل مسائل بهینه سازی بزرگ مقیاس پیوسته به کار 


الگوریتم بهینه سازی جامعه مورچه ها 
دردهه های اخیر استفاده از مفاهیم هوش جمعی 
(معصه‌عنااه‌اص1 ه/9۳) موجود در انا کی اجتماعی 
حشرات. به منظور ایجاد سیستمهای مصنوعی. توجه 
بسیاری از محققین را به خود جلب نموده است. از جمله 
این حشرات می توان به مورچه اشاره کرد. مورچه ها در 
یت رها کرد هل اما مس تقایل بسا اسساووران 
حس قوی بویایی خود کوتاهترین مسیر بین لانه و منبع 
غذا را تعیین نمایند. وقتی مورچه ای به دنبال غذا می 
گردد. در طول مسیر خود ماده بوداری به نام فرمون 
(۳۳۵۲۵۳۵۵86) به جای می گذارد که ساير مورجه‌هایی 
که در جستجوی غذا می باشند را به عبور از آن سیر 
تشویق می نماید. بر این اساس. 001070 و همکارانش 
در سال ۱۹۹۱ الگوریتم سیستم مورچه را معرفی نمودند 
[2]. به منظور حل مساله بهینه سازی با استفاده از 
الگوریتم بهینه سازی جامعه مورچه هاء در ابتداه می 
بایست گراف مناسب برای مساله تعربف نمود [3]. به 
عنوان نمونه گراف (),ا,)<0 را در نظر بگیرید. که 
[,......,یه ,)۲ مجموعه نقاط تصمیم می باشد 
که در این نقاط تصمیم می گیریم؛ (رب1]|-,آ مجموعه 
انتخاب های [ (10 ,....,1,2[) در هر یک از نقاط تصمیم 
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1( ...2 ,اح1) و 80 مجموعه هزینه های هر 
یک از انتخاب های مجموعه با میباشند. شکل (۱) نمونه 
ای از گراف تعریف شده برای مساله را نشان می‌دهد. در 
این شکل خط چین های عمودی نشاندهنده نقاط 
تصمیم دوایر توپر نشاندهنده گزینه های تصمیم و 
خطهای پررنگ مورب نشاندهنده یک جواب ساخته شده 


توسط مورچه فرضی می باشد. 


0 ...2 و1 


گزینه های تصمیم 


شکل ۱ گراف الگوریتم بهینه سازی جامعه مورچه ها 


بعد از تعریف گراف مناسب مساله. مورچه ها بر 
روی نقاط تصمیم مساله قرار داده شده و با انتخاب گزینه 
حرکت می کنند. انتخاب گزینه های تصمیم موجود در 
صورت تابع احتمالاتی ژر انبان امین کردد.: 


)۱( ]| )ره 


< (,6) رم 


ارجا اه 


از 
که در رابطه فوق. (4,) 0 احتمال اینکه مورچه1 
درتکرار ] و نقطه تصمیم ‏ گزینه تصمیم ژرا انتخاب 
نماید؛ (4) ریت مقدار فرمون موجود بر روی گزینه تصمیم 
[در نقطه تصمیم 1و تکرار ؛ ,1۲ اطلاعات کاوشی 
مربوط به گزینه تصمیم [در نقطه تصمیم ذو 0 و 8 
پارامترهایی هستند که نسبت وزنی بین فرمون و مقدار 
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بهره‌برداری بهینه از مخازن سدها با استفاده از مکانیزم .. 


کاوشی را کنترل می نمایند. 

با توجه به رابطه احتمالاتی فوق, مورچه ها در هر 
نقطه تصمیم. گزینه مورد نظر را انتخاب کرده و به نقطه 
تصمیم بعدی می روند. هنگامی‌که مورچه از کلیه نقاط 
نیقی وی کقلی یک مت انب یاف متسه فتاه 
است. این روند برای کلیه مورچه های یک جامعه در هر 
تکرار انجام می شود و قبل از شروع تکرار بعدی لازم 
است که فرمون مسیرها به روز رسانی (1/00216) شوند. 
رابطه به روز رسانی فرمون موجود بر روی گزینه های 
تصمیم در انتهای هر تکرار به شکل زیر تعریف می‌شود: 
۲( ()ر۵۲ +(1) ررعم ع (1+)) ره 


که در رابطه فوق؛ (1+1)ر,ت مقدار فرمون موجود بر 
روی گزینه تصمیم [در نقطه تصمیم ]در تکرار 1 
(1) رت مقدار فرمون موجود بر روی گزینه تصمیم [ در 
نقطه تصمیم 1 در تکرار 4 م ضریب تبخیر فرمون, که 
موجب کاهش میزان فرمون براساس مقدار ضریب 
تعریف شده می گردد. و (۸۳,)0 مقدار افزایش فرمون 
موجود بر روی گزینه تصمیم [ در نقطه تصمیم 1 در 
تکرار ؟ توسط مورجه > می باشد.. مقدار تغییرات فرمون 
( (۸۳:,)0) درمسائل حداقل سازی به صورت زیر تعریف 


اس نت 1 
مورچه >1 م انتخاب گردد 7 
1 
< () ۸۵۳ 
در غیر اینصورت 0 

۳( 
که در رابطه فوق (0) 18 مقدار تابع هدف به ازای جواب 
ساخته شده به وسیله مورحه 6 ام در تکرار ]و ۶ مقدار 
ثابتی است که ضریب با زگشت فرمون امیده می شود. در 
این مقاله» از بین الگوریتم های مختلف پیشنهادی جامعه 
مورچه هاء الگوریتم سیستم مورچه های بیشینه - کمینه 
انتخضاب شده و مسائل مورد نظر با بکارگیری این 
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الگوریتم مورد بررسی قرار گرفته اند؛ که در ادامه 


توضیحاتی مختصر درباره اين الگوریتم ارائه شده است. 


۳ 3 
الگوریتم سیستم مورچه های بيشینه - کمینه 
(۷۲۷۲۵۵6) 
الگوریتم سیستم مورچه‌های بیشینه - کمینه توسط 
۱ و 1005 [5] در سال 2۱۰ ارائه ح روت ایسن 
الگوریتم بر پایه اصول الگوریتم پایه اولیه پیشنهادی 
سیستم مورجچه ها استوار شده است؛ اما تغییراتی به منظور 
بهبود روند همگرائی در الگوریتم سیستم مورچه ها 
صورت گرفته استت: در الگوریتم سیستم مورچه های 
بیشینه - کمینه با محدود کردن مقادیر فرمون به یک 
محدوده نعیین شده تس از همگرایی بیج آن 
اجتناب می شود. مقادیر محدوده بالائی و پائینی فرمون 
(فرمون حداکثر و حداقل) در تکرار ] توسط روابط زیر 


ماب فی آر 03: 

4 1 
( 2-0( هو ۲ 
طلست 


هل (1- ,0 )1( زو 


که در روابط فوق؛ (() *)2 مقدار تابع هدف در بهترین 
مسیر طی شده توسط مورچه ها تا تکرار ! م ضریب 
تبخیر فرمون؛ 1 ضریب بازگشت فرمود؛ ,وه 0 
ضریبی است که هر چه مقدار ان کمتر باشد مرزهای 
محدوده فرمون ریزی به هم نزدیک شلده و به تبع آن 
احتمال انتخاب گزینه های تصمیم تقریبا نزدیک به هم 
می شود و ,,, 0 میانگین تعداد گزینه های تصمیم 
موجود در نقاط تصمیم می باشد. به این ترتیب در انتهای 
هر تکرار مقدار فرمون مسیر ها کنترل می گردد تا از 
مقدار (1) ,,,> تجاوز نکنند و از مقدار (۶) مزیرت 
کمتر نگردند. همچنین این الگوریتم از رابطه زیر به 
منظور به روز رسانی فرمون مسیرها استفاده می کند: 
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1۹ 


0 (۱) :۸۲۳ +()) عم <(1+)) ره 

در رابطه فوق, (۶) :۸۳ فرمون اضافه شده به هر 
یک از جزء مسیرهای (گزینه تصمیم [ در نقطه تصمیم 1) 
بهترین مسیر طی شده توسط مورچه ها (بهترین مسیر 
توسط مورجه 1 طی شده و جواب ,5 را ساخته است) 


در تکرار ] ام می باشد که از رابطه زیر محاسبه می شود: 


۱ 0 


که در رابطه فوق. (2) ,1 تابعی است که مقدار آن به 
ازای ۸ 2 برابر با یک و در غیر این صورت برابر با 
صفر می باشد و بقیه پارامترهای رابطه پیش از این 


تعریف شده اند. 


مکانیزم تظریف تطبیقی قطعی 
از جمله ویژگی‌های الگوریتم پایه اولیه پیشنهادی جامعه 
مورچه‌ها می‌توان به حرکت مورچه های مصنوعی در 
فضای گسسته اشاره نمود. ماهیت این الگوریتم پایه اولیه 
جستجو در فضای جستجوی گسسته مساله را به این 
رمک ام اماب ستاری از سای واقس دای 
اطراف ما پیوسته می باشند و بنابراین به منظور حل این 
مسائل واقعی پیوسته با استفاده از این الگوریتم لازم 
شین که سرام کشوم تسا له تسا کسستاه متا رمع 
محدوده های متغیرهای تصمیم مساله. تشکیل شود. 
چگونگی گسسته سازی متغیرهای تصمیم مساله بر روی 
کشت ربیاو کیت هزات مساله تان.می کذاره 
به گونه ای که با گسسته سازی درشت متغیرهای تصمیم 
مساله ممکن است نقاط بهینه در فضای جستجو مساله 
قرار نگیرند و گسسته سازی ریز موجب بزرگ شدن ابعاد 
مساله و افزایش تلاش محاسباتی و بعضا کاهش کیفست 
جواب می گردد. بنابراین. افزایش قابلیت‌های الگوریتم 
بهینه سازی جامعه مورچه ها به منظور نزدیکتر شدن 


تشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


۷/۰ 


جوابها به جوابهای بهینه» در حل مسائل بهینه سازی 
پیوسته. یکی ازچالشهای این الگوریتم محسوب می شود. 
بدین منظور, افشار [14] مکانیزم مناسبی به نام مکانیزم 
تظریف تطبیقی قطعی ([1۸) به منظور رسیدن به جواب 
هایی نزدیک به جواب های بهینه مطلق مساله ارائه نمود. 
این مکانیزم. مکانیزمی تکراری است که در آن فضای 
جستجو مساله در هر تکرار به اطراف جواب بهینه ایجاد 
شده در جستجوی قبلی محدود می شود. به این ترتیب 
با محدود کردن فضای جستجو مساله. فاصله گزینه همای 
تصمیم در هریک از نقاط تصمیم از یکدیگر کمتر میگردد 
و بنابراین پس از چند بار اعمال این مکانیزم. جواب‌های 
حاصله به جواب‌های بهینه مساله واقعی نزدیک‌تر 
می‌گردند. این روش جایگزین بسیار خوبی به جای 
گسسته سازی های زیاد متغیرهای تصمیم مساله و 
حستجوی در فضای جستجوی بزرگ ساخته شده مساله 
می باشد. همچنین. این روش نیاز به هزینه عملیاتی 
کمتری دارد و جواب‌های دقیق‌تری نیز با بکارگیری آنه» 
برای مساله حاصل می‌شود. در این مقاله از دو روش 
برای محدود کردن فضای جستجوی مساله و اعمال این 
مکانیزم استفاده شده؛ که در ادامه توضیحاتی درباره اعمال 
هریک از آنها ارائه شده است. 

روش اول: در این روش محدوده همای فضای 
جستجوی جدید مساله شامل تعداد مشخصی از 
گزینه‌های تصمیم تعریف شده اولیه در همسایگی جواب 
بهینه یافته شده در جستجوی قبلی می باشد. در شکل 
(۲) فضای جستجوی جدید ساخته شده برای مساله با 
اعمال روش اول مکانیزم. که در ساختن فضای جستجوی 
جدید مساله در هر یک از نقاط تصمیم از یک گزینه 
تصمیم در همسایگی جواب بهینه یافته شده قبلی استفاده 
شده نشان داده شده است. در این شکل ابتدا جواب 
بهینه یافته شده درجستجوی قبلی نشان داده شده است 
(شکل ۲-لف). سپس یک گزینه تصمیم در همسایگی 
جواب بهینه یافته شده قبلی. به منظور تعیین محدوده 
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بهره‌برداری پهینه از محازل سدها با استنفاده از مکانیزم ی 


های جدید متغیرهای تصمیم و ساختن فضای جستجوی 
جدید مساله» انتخاب گردیده است (شکل ۲- ب». در 
نهایت در شکل (۲-ج) فضای جستجوی جدید مساله 
که با گسسته سازی محدوده جدید بدست آمده برای هر 
متغیر تصمیم و تعیین گزینه های موجود برای آن ساخته 
می شود. نشان داده شده است. در این شکل ها خحط 
چین‌های عمودی نشان‌دهنده نقاط تصمیم دوایر توپر 
نشان‌دهنده گزینه های تصمیم. خطهای پررنگ مورب 
نشان‌دهنده جواب ساخته شده توسط مورچه فرضی و 
خط چین های مورب نشان دهنده محدوده های جدید 
تعریف شده برای هر متغیر تصمیم می باشد. 

۲ روش دوم: در این روش مرزهای جدید محدوده های 
متغیر تصمیم مساله به منظور ساختن فضای جستجوی 
جدید مساله. با تعبین یک بازه اختیاری شامل جواب 
ی نا ی ترس مضه و 
آن دامنه بازه از طریق کاربر تعیین می شود. در شکل (۳) 
چگونگی انجام این روش در اعمال مکانیزم گسسته 
سازی مجدد متغیرهای تصمیم مساله. به منظور ساختن 
گراف جدید مساله نشان داده شده است. در ایین شکل 
ابتدا جواب بهینه یافته شده در جحستجوی قبلی نشان داده 
فبت اس (قکل اتف )سس ینک سازه ارم 
شامل جواب بهینه بدست آمده از مرحله قبل به منظور 
ساختن فضای جستجوی جدید مساله انتخاب شده است 
(سکل ۳- ب). در نهایت در شکل (۳-ج) فضای 
جستجوی جدید مساله. که با گسسته سازی محدوده 
جدید بدست آمده برای هر متغیر تصمیم و تعیین 
گزینه‌های موجود برای آن ساخته می شود نشان داده 
شده است. در این شکل ها خط چین های عمودی 
نشان‌دهنده نقاط تصمیم, دوایر توپر نشاندهنده گزینه‌های 
تصمیم. خط‌های پررنگ مورب نشاندهنده یک جواب 
ساخته شده توسط مورچه فرضی و خط چین های 
مورب نشان‌دهنده محدوده های جدید تعریف شده برای 
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۷۱ 


٩ ۵ : ۳ ۰ ۰ ۲ 0‏ 
۱ ۰ ۱ ۱ ۱ ۱ ۱ ۱ ۱ ۱ 
۱ ۰ ۰ ۵ . ۰ ۰ ۰ ۰ 8 
ج ۵۵ 2 ۲ 9 ه و ها 
ک ا ره ۱ 0 
‌ ۱ ‌ ۰ .۹ ۱ 21 
۰ ۱ ۰ ۱ 9 
۰ ۰ ه‌ ۰ رگررینه ۰ 

‌ 0 ح[ 

شکل ۲ نمایش شماتیک روش اول اعمال مکانیزم *۸1 به منظور ساختن گراف جدید مساله 
ها و 0 ۳ ۶۰ 9 ۲۰ 
4 ۱ ِِِ ۳۹ [ ۰ ۰ ۱ 
۱ ه: ۱ 1 و را ۱ ۱ ۳4 ۷ 

۱ ۰ . ۰ / ه ٩‏ .8 
و . و ح ی 
! و ۲ ۱ 0 
و ۱ ۳ ۲ ۳ ٩‏ یه ٩‏ 

,0 0 ع[ 


شکل ۳ نمایش شماتیک روش دوم اعمال مکانیزم ۸8 به منظور ساختن گراف جدید مساله 


مسائل نمونه 
به منظور بررسی قابلیتهای مکانیزم های پیشنهادی. مسائل 
نمونه‌ای با اعمال روش‌های اول و دوم مکانیزم ۲۸۶ و 
بکارگیری الگوریتم سیستم مورچه های بیشینه - کمینه 
حل یله استت. لذا در این تحقیق. در ابتداء دو تابع 
ریاضی بیچیده 7 و 2601 مورد بررسی قرار گرفته 
است و در ادامه مساله بهره برداری بهینه ساده و برقابی از 
مخزن سد دز بعنوان یک مساله بهینه سازی واقعی در 


حوزه مهندسی اب. مورد بررسی قرار گرفته است. 


تایع ۸۰/6 . به عنوان اولین مساله نمونه ریاضی. یافتن 
حداقل مقدار تابع ریاضی ۸۵116۷ مورد نظر می باشد. 
تابع 7 یک تابع پیوسته با توابع غیر خطی نمایی و 
کسینوسی می‌باشد. شکل کلی تابع مذکور به صورت 
زیرمی باشد [16]: 
| و 02 اس ۵ +20 (۲0 
2 


رماع کک > 5- | ( مرت اجه 
21 1 


(۸) 
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این تابع دارای جندین حداقل موضعی است که از 
هر طرف به سمت یک حداقل مطلق در حرکت 
ازای منظور کردن مقدار صفر برای متغیرهای تصمیم 
حاصل می شود. این تابع در حالت کلی دارای 7 بعد می 
باشد که در این تحقیق شکل (۱۰) بعدی این تابع برای 
به دلیل ماهیت تصادفی آن‌ها دارای پارامترهای آزادی 
هستند که با تغییر مقدار این پارامتر ها؛ عملکرد الگوریتم 
همگرایی و مرغوبیت جواب‌های حاصل از به کارگیری 
آزتهتا ی هراشا که یراق شم مسا لها با رف 
سعی و خطاء مقادیرمطلوب پارامترهای الگوریتم ۱/۵ 
به صورت مقادیر جدول (۱) ارائه شده اشنتت: در این 
مساله» با توجه به شکل تابع امکان ارائه تعریفی 
مناسب به منظور محاسبه مقادیر کاوشی وجود نداشته و 
بنابراین مقدار ضریب تاثیر آن صفر (8<0) منظور شده 


است. 
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جدول ۱ مقادیر مطلوب پارامترهای الگوریتم ۱۸۸5 


تعداد مورجه‌ها 


0 0 9۷ امه ۳ 
۷ ۱ ۰ ۰۵ ۰/3 


که در بازه [ ۵0,۵-] تعریف شده است. با استفاده از تقسیم 
بندی ۰۱۰ ۰ و ۰ تایی گسسته سازی شده است. نتایج 
حاصل از ۱۰ بار اجرای مدل. با استفاده از روش اول 
اعمال مکانیزم ۸8 و بدون استفاده از آن در جدول 
(۲) ارائه شده است که در آن یک بار از یک و بار دیگر 
از دو گزینه تصمیم در همسایگی جواب بهینه یافته شده 
قبلی به منظور تعریف فضای جستجوی جدید مساله 
استفاده شده است. حداکثر تعداد تکرارها در تمامی 


اجراها برابر ۱۰۰۰ در نظر گرفته شده است. همانگونه که 


بهره‌برداری بهینه از مخازن سدها يا استفاده از مکانیزم .. 


از مقایسه نتایج جدول (۲) مشاهده می شود بدون 
انتتفاته از مکانیزم لفط در سیم بتدی رفس 
متغیرهای تصمیم مساله دقت جوابهای یافته شده بسیار 
کم می باشد و تقسیم بندی ریز متغیرهای تصمیم مساله 
نیز نتوانسته است به اندازه اعمال مکانیزم 9۸8 در یافتن 
جواب مناسب موثر واقع شود. همچنین در حالت بدون 
استفاده از مکانیزم 0۸ الگوریتم دچار همگرایی زود 
شش له یز رافونقن تدش سانشان فتاییی فتاه 
تکرارها و يا مورچه ها) نیز نمی تواند منجر به یافتن 
جواب‌های مناسبتر گردد. در حالی که اعمال مکانیزم 
۶ موجب ایجاد شک در مکانیزم جستجو و فرار از 
گیر افتادن جواب‌ها در نقطه بهینه محلی و حرکت به 
سمت جواب های مناسبتر (نزدیک به بهینه مطلق) 


می‌شود. 


جدول ۲ نتایج حداقل سازی تابع ۸۰۸6۷ بدون استفاده و با استفاده از روش اول اعمال مکانیزم ۸ 
بدون استفاده از مکانیزم ۱3۹4 با استفاده از روش ول مکانیزم ۲۸۶ 
منت یو شش | ۲۰ شیم دی خ 
یک گزینه تصمیم | دو گزینه تصمیم 
۱ 2۱2-۳۰۷۹۵ ۱/۸۳۵۹۹۵ ۳۳۳۲۳ ۹۳« «/,«",۱,, ۱/6 
۲ ۳۰۷۵:-2 ۱۵ ۳۳۳۲۳ ۱۳۰۹۹ ۰۹ 
۳ 32326 ۱/۹ ۰۳۳۳۳ ۱۳2 ۰۱۳۵5۹ 
3 3236 ۱/۹ ۰۳۳۳۳ ۱۳۸2۸ ۳/۳۲ 
۵ 2-2-۳۷۹۵ و ۳۳۳۳ ۱۱۳۳" ۳/۲ 
1 2۱(۳۰۷۵ ۱/۸۹۵ ۳۳۳۳ 2 ۰/۱۳۵ 
۷ 2-2-۳۷۹۵ ۱۱/۸۳۹۵ ۳۳۳۳ 2 ۱۹ 
۸ 2(۳۷۵ ۱۳۸۹۵ ۳۳۳۳ ۱/۳/۳۹ 32 
٩۹‏ 2۱(۳۰۷۵ ۱۸۹۵ ۳۳۳۳ ۳۸۱ ۳/۵۸ 
۰ 2۱(۳۷۵ ۱/۸۹۵ ۳۳۳۲۳ ۹۱*:۹/:۵/ ۰/۱۰۹۹ 
ميانگین 12-1۳۷۵ 0 ۳۳۳۳ ۱۹۳۳۷ ۱/۳۹۵ 
بهترین جواب 21:۳۵ ی ۳۳۳۳ ۰۹/۹۳+«,«-1-۰ ۹ 
بدترین جواب 2۱(۳۰۷۵ ۱۸۹۵ ۳۳۳۳ ۳۸۱۹ 2 
انحراف معیار ۰ ۰ ۰ 1۳ ۱۹۹2۷ 
ضریب تغییرات ۰ ۰ ۰ ۱۴ ۰۹۸۳ 


نشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


سال بیست و سوم شماره یک. ۱۳۹۰ 


محمدهادی افشا ر- ابراهیم رضایی سنگدهی- رامتین معینی 


۴ که تابعی غیر خطی و نمایی است. مدنظر قرار 
)+ )+2۳ ری؟ 
9, و 2 ۱[ 1۳ 1> ۶ > 1- 
)۹( 
(( رم رو > 1 ۵ +۱ م1 )وه +1 - و 
1 (ر رک )9۵ )۲ یک 
03 
- (( )۵ +1)/ ر«<1-) مه +1 و 
( رم ((ر )۵ +1 / )۳ 
۱۱2( 
- ((رج <1-)00 +1)/ س ۷-) مه +1 9 
( رم - (ار«- ری )00 +1)/ )۳ 
ره 
هه ون > <1 ۲ +۱۱ ت «-) مه +۱ - و 
( رو - ((را- <1-)90۳0 +70 و 


(۱۳ 


مساله مورد نظراز نوع حداکثر سازی و دارای ٩‏ 
متغیر تصمیم است. برای این مساله نیز مقادیرمطلوب 
پارامترهای الگوریتم ۵۹ به صورت مقادیر جدول 
(۱) ارائه شده است. به منظور حل مساله» متفیرهای 
تصمیم مساله» که در بازه [۱0۱-] تعریف شده است. با 
استفاده از تقسیم بندی ۱۰ تائی (تقسیم بندی درشت) و 
۰ تایی (تقسیم بندی ریز) گسسته سازی شده است. 
نتایج حاصل از ۱۰ بار اجرای مدل با استفاده از روش 
اول اعمال مک‌انیزم 9۸۴ و بدون استفاده از آن» در 
تور :(۳) اراکه: شده: ات که:در ان شنک سار از سکاو 
بار دیگر از دو گزینه تصمیم در همسایگی جواب بهینه 
یافته شده قبلی به منظور تعریف فضای جستجوی جدید 
مساله اتتفاوم فته سک سا | کر مداد نکر ارت ور 
تمامی اجراها برابر ۱۰۰۰ در نظر گرفشه شسده امست: 
همانگونه که از مقایسه نتایج جدول (۳) مشاهده می‌شود. 
در حالت استفاده از روش اول اعمال مکانیزم ۲۸ 
جواب های مطلوبتری نسبت به حالت بدون استفاده از 
مکانیزم 1۸1 در هر دو حالت تقسیم بندی درشت و 
ریز متغیرهای تصمیم مساله به دست آمده است. 


جدول ۳ نتایج حداقل سازی تابع 608 بدون استفاده و با استفاده از روش اول اعمال مکانیزم 0۸1 


سال بیست و سوم شماره یک» ۱۳۹۰ 


تشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


بدون استفاده از مکانیزم ۳۸۶ با استفاده از روش اول مکانیزم ۲۸۶ 
وت هورق ای و 
یک گزینه تصمیم | دو گزینه تصمیم 
۱ ۲ ۱/۷ ۱/۹۲ ۱/۹۲ 
۲ ۲ ۱/۹۸ ۳/۹۲ ۳/۹۲ 
۳ ۱/۹۸ ۱/۹۹ ۲/۳۹ 1/۸ 
3 ۲ ۱/۹۹ ۲ ۲ 
۵ ۲ ۱/۹۹ ۳/۹۲ ۳/۹۲ 
1 ۲ ۱/۹۹ ۳/۹۲ ۳/۹۱ 
۷ ۲ ۱/۹۹ ۱/۳۵ ۱/1۹ 
۸ ۲ ۱/۹۲ ۱/۳۷ ۱/۷ 
۹ ۱/۹۹ ۱/۹۹ ۳/۹۲ ۳/۹۲ 
۱۰ ۲ ۲ ۱/۹۲ ۱/۹۲ 
میانگین ۲ ۱/۹۸ ۱/۹ ۳/۳۷ 
بهترین جواب ۲ ۲ ۳/۹۲ ۱/۹۲ 
بدترین جواب ۱/۹۸ ۱/۹۲ ۲ ۲ 
انحراف معیار ۰«۰«۹91۸۲,+* ۹۳۳+ ۱۷-+(-2,+* 1 
ضریب تغییرات ۰«.«(۵, ۰(۹",«,* 2 (2 


۷ 


مساله بهره برداری بهینه از مخزن سد. به منظور نشان 
دادن قابلیت‌های مکانیزم 9۸ در حل مسایل حوزه 
مهندسی آب؛ دو مساله واقعی تک هدفه بهره برداری 
بهینه ساده و برقابی از مخزن سد مورد بررسی فرار 
گرفت. اين مسائل از جمله مسائل بهینه سازی مقید می 
باشد, که به منظور تعریف مساله. لازم است که تابع 
هدف و قیود مساله تعریف شوند. تابع هدف در مساله 
بهره برداری از مخزن می تواند به صورت‌های مختلف 
تعریف گردد. تابع هدف می تواند حداکثر سازی درآمد» 
حداقل سازی هزینه. حداکثر سازی انرژی تولیدی؛ 
حداقل سازی میزان کمبود اعمال شده به یک نیاز 
ی ور کرک در دراب یو ور 
هدف حداقل سازی میزان کمبود اعمال شده نسبت به 
یک نیاز مشخص (در مساله بهره برداری ساده) و حداقل 
سازی کمبود توان تولیدی نسبت به ظرفیت نصب 
نیروگاه (در مساله بهره برداری برقابی) مدنظر قرارگرفته 
است. 

به عنوان مسائل نمونه واقعی. از اطلاعات سد دز 
استفاده شده است. سد دز بلندترین سد ایران و در زمان 
ساخت نیز مرتفع ترین سد مخزنی در خاورمیانه بوده و 
بر روی رودخانه دز در استان خوزستان. جنوب غربی 
تایه کید ایتک سا سس وی ۲۵ عرش شا 
شهرستان دزفول قرار دارد. رودخانه دز که از ارتفاعات 
غربی زاگرس (کوههای بختیاری) سرچشمه می گیرد از 
نظر میزان آبدهی دومین رودخانه ایران محسوب می شود 
و در 1۵ کیلومتری شمال اهواز به رودخانه کارون 
فص اد نی ی 
ارتفاع آن از پی ۲۰۳ متر و از کف رودخانه ۱۹۰ متر 
است. عرض بدنه در پی ۲۷ متر و در تاج ۶/۵ مت طول 
تاج ۲۱۳ متر و رقوم تاج سد ۶ متر از سطح دریا 
می‌باشد. سطح دریاچه در این رقوم به 1۵ کیلومتر مربع 
می‌رسد. 

در این مساله. حجم حداقل و حداکثر مخزن به 
ترتیب معادل ۸۳۰ و۳۳۶۰ میلیون متر مکعب و حجم 


نشریه مهندسی عمرال دانشگاه فردوسی مشعد 


بهره‌برداری پهینه از محازل سدها با استفاده از مکانیزم 


اولیه مخزن ۱۳۶۰ میلیون متر مکعب می باشد. از بین ۰؟ 
سال آمان ۵ سال ابتدایی (۲۰ دوره ماهیانه) به منظور 
مدلسازی انتخاب گردید. مقادیر نیاز و ورودی های 
ماهیانه مشخص می باشد. به گونه ای که جمع ورودی 
ماهیانه به مخزن در این ۰" ماه معادل ۲۵۱۵ میلیون متر 
مکعب و جمع نیاز ماهیانه معادل ۲۹۵۰۲ میلیون متر 
مکعب میباشد. مقدار جریان ورودی به مخزن سد و 
مقادیر نیاز ماهانه در این ۵ سال آمار در شکل (4) ارائه 
شده است. بدین ترتیب مقدار مجموع باز از مجموع 
آورد در این ۵ سال به اندازه ۲۹۸۲ میلیون متر مکعب 
بیشتر میباشد. علاوه بر آن» در این مساله. جهت تعین 
مقدار بار آبی موثر بر توربينها. تراز پایاب نیروگاه سد دز 
معادل ۲ متر از سطح دریا منظور شده است. این سد 
۰ هکتار از اراضی پایین دست را آبیاری می کند و 
نیروگاه آن از ۸ واحد ۸۰/۸ مگاواتی تشکیل شده است 
که زمان کارکرد آن در طول روز حدود ۱۰ ساعت 
می‌باشد. به این ترتیب در محاسبات ضریب کارکرد (۳۳) 
معادل ۰/۶۱۷ منظور شده است. ظرفیت نصب پیروگاه 
معادل 1۵۰ مگاوات و بازده آن معادل ۸٩۰‏ منظور گردیده 


تیت: 


نیاز سس 5 1800 


جریان ورودی - سب - ۱ 1600 
14001 
۱ 1 1200 


1000 1 ۱ [ 


نیاز و جریان ورودی ماهانه (/0/01) 


0 55 50 45 40 35 30 ۰.25 ۰ 20 15 10 5 0 
دوره زمانی 


شکل ۶ مقادیر جریان ورودی و نیاز ماهانه در دوره زمانی ۵ ساله 


مساله بهره برداری ساده از مخزن ۳ ۱۷۶ در این مساله 
تابع هدف به صورت حداقل سازی میزان کمبود اعمال 


شده نسبت به یک نیاز مشخص تعریف شده است: 


13 ۱ [ی/(-,0)] 02 مت صنصر 


احا 


سال بیست و سوم شماره یک. ۱۳۹۰ 


محمدهادی افشار- ابراهیم رضایی سنگدهی- رامتین معینی 


که در رابطه فوق. , ٩‏ مقدار نیاز در دوره زمانی ] 
(ماهیانه)؛ ,7 مقدار آب رها شده از مخزن در دوره زمانی 
او ره حداکثر نیاز در تمام دوره‌های زمانی (7۲) 
می‌باشند. 

اولین گام به منظور حل مساله بهینه سازی با کمک 
الگوریتم بهینه سازی جامعه مورچه هاء تعریف گراف 
مساله می باشد. بدین منظور لازم است که نقاط تصمیم. 
متغیرهای تصمیم و گزینه های متغیرهای تصمیم مشخص 
شوند. در این مساله, دوره های زمانی بهره برداری از 
مخزن به عنوان نقاط تصمیم و مقدار آب رها شده از 
مخزن به عنوان متغیر تصمیم منظور شده است. گزینه 
های تصمیم مساله برابر با گسسته سازی های صورت 
گرفته برای مقادیر آب رها شده از مخزن سد در دوره 
های زمانی می باشد. شکل (۵) گراف مساله بهره برداری 
از مخزن. در حالتی که مقدار آب رها شده از مخزن 
بعنوان متغیر تصمیم منظور شده است را نشان می دهد. 

اکثر مسایل بهینه سازی مقید می باشند. بنابراین 
لازم است علاوه بر تابع هدف قیود مساله نیز تعریف 
شوند. قیود موجب می گردند تا برخی از جواب‌ها غیر 
شدنی باشند و بنابراین لازم است با اعمال مک‌انیزمی 
مورچه ها را به جستجو در فضای شدنی ترغیب نمود. 
قیود مساله بهره برداری ساده از مخزن به صورت رابطه 
پیوستگی, حداکثر و حداقل مقدار آب رها شده از مخزن 
و حجم مخزن تعریف می گردند. به طور کلی می توان 
قیود مساله را به شکل زیر تعریف نمود: 


۷۵ 


۱۷ رگ اک 5 


که در روابط فوق. 5 مقدار حجم مخزن در دوره 
زمانی ؛ ,۹ مقدار حجم مخزن در دوره زمانی ۱+ 
# , 4 مقدار اب ورودی به مخزن در دوره زمانی 4 ,۲ 
مقدار اب رها شده از مخزن در دوره زمانی 5 هت 
حداقل مقدار اب رها شده از مخزن؛ یی * حداکتر مقسدار 
اب رها شده از مخزن؛ 9 حداقل مقدار حجم مخزن 
ق نس 5 حداکثر مقدار حجم مخزن می باشد. در این مساله 
از روش تابع جریمه برای اعمال قبود استفاده شده است. 
در این روش تابع هدف مساله بصورت زیر اصلاح 
می‌شود: 
۱ ۱ ۷ 
هه رت 4 0-2 
احا احا 


که در آن 
نگ ک 5 1 زگ / (5- )٩«‏ >« < ۳00 
٩‏ < 5 1 مگ 2066 ۰ ۳۵0 
)۱۹( 


که در روابط فوق. موم و "موم به ترتیب مقدار جریمه 
در نظر گرفته شده برای قید حداقل حجم و حداکثر 
حجم مخزن در دوره زمانی ]و 6 ضریب ثابتی» که در 
این مساله برابر با ۱/۷ در نظر گرفته شده است. می باشد 
و سایر پارامترها پیش از این تعریف. شده اند با ترخنه 
به تابع هدف مساله مقدار اطلاعات کاوشی برای هر یک 
از جزء مسیرها (گزینه های تصمیم) بصورت رابطه زیر 
بر حسب مقدار تابع هدف برای هر یک از جزء مسیرها؛ 


تعر یف شو د: 
(۱۵( 25 رب5 7 ٍِِ«« 2 1 
:1۱ 

۱10( مگ گِ 1 ‌ِ م هط [ »»۵/( 7 7 

۶۶ ۰ ۰ ۶ 

ا ا ۱ 

‌ ‌ ‌ جل) مب مد 

هب هب اه 3 

۰ ۰ ۰ وخ 


نقاط تصمیم ( دوره های زمانی بهره برداری) 


- 


شکل ۵ گراف مساله بهره برداری از مخزن 


سال ب بیست و سوم شماره یک» ۱۳۹۰ 
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۷۹ 


که در رابطه فوق. ره ژامين گزینه تصمیم (مقدار آب 
رها شده از مخزن) در نقطه تصمیم 1 (دوره زمانی 1 4 
مقدار نیاز در نقطه تصمیم 1 (دوره زمانی ) و ,مه برابر 
با حداکثر مقدار نیاز در کل دوره های زمانی می باشد. به 
منظور بررسی نحوه عملکرد الگوریتم جامعه مورچه ها و 
مکانیزم پیشنهادی (مکانیزم؟۸1). در مساله بهره برداری 
ساده از مخزن سد. از اطلاعات سد دز استفاده شده است. 
در این تحقیق, به منظور تعریف گراف مساله. محدوده 
مجاز مقدار آب رها شده از مخزن [۱۰۰۰۰۰] می باشد. 
جدول (4) مقادیر مطلوب پارامترهای بکار گرفته شده در 
الگوریتم 10۸۹ به منظور حل مساله بهره برداری ساده 
از مخزن سد دز را نشان می دهد. 


جدول ۶ مقادیر مطلوب پارامترهای بکار گرفته شده در الگوریتم 
۱/۸۵۹ 


تعداد مورچه ها 0 8 0 امه ۲ 
۱۰۰ ۱ ۰/۳ ۰۲ ۰/۹۵ 


حدول (۵) نشان دهنده نتایج ستفاده از روش اول 
اعمال مکانیزم 9۸1۴ در حل اين مساله می باشد که در 


بهره‌برداری بهینه از مخازن سدها يا استفاده از مکانیزم .. 


آن از تعداد گزینه های تصمیم متفاوت. به منظور تعریف 
فضای جستجوی جدید مساله» بهره گرفته شده است. در 
۰ درنظر گرفته شده است. همانگونه که از مقایسه 
نتایج حدول 6 مشاهده می‌ شود استفاده از دو کربفه 
تصمیم در همسایگی جواب ساخته شده قبلی, به منظور 
تعریف فضای جستجوی جدید مساله در روش اول 
اعمال مکانیزم ۴ به نتایج بهتری منتهی شده آسستتا: 
همچنین. در جدول (1) نتایج روش دوم اعمال مکانیزم 
وبا یره کر از بتاژه های اف اراته شبن 
است. مقایسه نتایج نشاندهنده انست که با انتخاب بازه 
۸ نسبت به ساير بازه ها؛ جواب مطلوب‌تری برای این 
مساله حاصل عابتا همچنین نتایج بهترین جوابهای 
حاصل از استفاده از روش اول و دوم اعمال مک‌انیزم 
۶ با جوابهای حاصل از بکارگیری الگوریتم ۸۸۱۸5 
بدون اعمال مکانیزم ۱۳3۹۹4 در جدول 02 ارائه شده 
است. مقایسه نتایج نشان می‌دهد که استفاده از مکانیزم 
تاثیرقابل ملاحظه ای در بهبود جواب‌ها داشته 


الشتن 


جدول ۵: مقایسه نتایج روش اول اعمال مکانیزم ۶ به ازای تعداد گزینه های تصمیم متفاوت 


تعداد اجرا 
۱ 
۱ 1 
۲ 0 
۳ ۱۹۹ 
3 1۹-۳۹ 
۵ ۱۳۱۷ 
1 ۱۹۳۸ 
۷ ۱۹ 
۸ ما 
۹ ۳ 
۷۰ ۱۷ 
میانگین ۱۹۸۶ 
بهترین جواب | ۰۸۷۹۱۹۰۲ 
بدترین جواب | ۲/۳۱۲۱۷ 
انحراف معیار ٩۸‏ 


ضریب تغییرات | ۰/۳۷۳۵۹۱ 


نشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


با استفاده از روش اول مکانیزم ۲۸۶ 


تعداد گزینه های تصمیم 
۲ ۳ 3 
۸۶ | ۰/۷۲۹۶ 1 ۰/۸۲:۳۰۳ 
۱ ۰/۹۸۰ 
۷ | ۰/۹۱۱۰۱ 1 ۰/۲۶۱۹ 
۸۲ | ۰/۸۲ | ۰/۸۸۹۸ 
۵ | ۱/۹۰۷۸ [ ۱۸۰۹۱۹ 
۹ ۱/۹۰۷۸ [ ۱/۲۵۸۷ 
۳ | ۱/۶۵۱۵ [ ۱/۰۷۲۷ 
۱ ۱ ۰/۹۱۹۸۷۲ 1 ۰/۹۲۹۲۰۲ 
۵۲۲۸ | ۰/۹۹۰۳ 
۷ | ۱/۰۳۰۷ [ ۱۰۳۶۰۷ 
۵ | :۱/۰۵۲۲ [ ۱۸۰3۹۵۵ 
۸۶ | ۰/۷۲۹۶ ]1 ۰/۸۲:۳۰۳ 
۵ | ۱/۹۰۷۸ [ ۱/۸۹۱۹ 
۱ ۱ ۰۳۰۹۲۹ 1 ۰/۲۸۵۸۹۸ 
۵ | ۰۰۹۵۳۵ ۰/۲۹۷/۳۰۷ 


سال بیست و سوم شماره یک. ۱۳۹۰ 
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جدول 1 مقایسه نتایج روش دوم اعمال مکانیزم به ازای طول بازه های متفاوت 


با استفاده از روش اول مکانیزم ۸ 
تعداد اجرا طول بازه 
۰/۳ ۵/ ۰/۷ ۰/۸ ۰/۹ 
۱ ۱۱/۰۹/۳۷ 2 ۰*۵۸ ۲ ۳۵ ,+ 
۲ ۱۱۰۷۲ ۱۱/۱/۰۳ ۱۳2 ۰,۱2۳ ۱۱/۰۷ 
۳ ۱۳۸۰ ۱۳2 ۱۱-۸۸ ۱2 ۱۳۸ 
1 ۱۱۳۳۹۵ ۵ ۰/۶۰۸۵ ۱/۰۰۵۹ ۰/۹۹۹۵ 
۵ ۱2۰3 ۱۳۸2 ۱۱۳۸ ۱۱,۸۶۹ ۱۵۵۵۳۸ 
1 ۱۳۶" ۱۳۸ ۱۸۵۵ ۱۱۳۱۷ ۱۱۶۸ 
۷ ۱۵ ۱۱۸۳/۳۵ ۱۱/۱۳۸ ۱۱۸۳۰۸۷ ۱۳2 
۸ ۱/۹۹۲۱ ۱/۳/۳۹ ۱/۰۳۸ ۱/۷۵ ۱/۰۹۲ 
۹ ۱-۸۰۱ 2( ۹1 ۳۲" ۰۱۸۳۷۲ 
۷۲۰ ۷ / ۷ ۷ ۱-۱ ۱۳/۷ 
میانگین ۸ ۳ ۱/۲۳۹۹ | ۱/۰۹۱۳ 1 ۱/۱۲۵۰ ۱۳۳/6۹ 
بهترین جواب ۱۳۸۵ ۳۹ ,+ ۰*۵۸ ۲ ۱( ۸" 
بدترین جواب | ۲/۵۶۲۷۹ | ۱۷۵۰۹۵ | ۱/۳۵۹۸ | ۱/۸۸۵۵ | ۲/۲۵۱۳۸ 
انحراف معیار | ۰/۸۲۲۳۶ ]| ۰۱/۲۲۷۵۲۱ | ۰/۲۳۰۵۱۲ | ۰/۳۰۱۳ | ۰/۰۵۵۲۳ 
ضریب تغییرات | :۰/۳۷۲۲ | ۰/۱۸۳۵۳۲ | ۰/۲۱۰۲۵۸ | ۱/۲۱۷۹۷۳ ۱۱ 
جدول ۷ نتایج مساله بهره برداری ساده مخزن بدون استفاده و با استفاده از روش اول و دوم عمال مکانیزم 1-۸ 
۱ بدون استفاده از با استفاده از مکانیزم ۸ 
تعداد اجرا 
مکانیزم ۳۸ | روش اول روش دوم 
۱ ۳۰۵۹۵(-,* ۰*۱۵ ۲ ,+" 
۳۲ ۱,(۳۵۹۵,* ۳۷" ۰,۱۳۳ 
۳ ۱۳۸ ۱,۵۳۷,* ۱۱/۰۹۸ 
3 ۱۳/5۷ ۱-۸۵ ۱/6۹ 
۵ 2 ۱۸,۵ ۱/۸۳/۸۹۶۹ 
1 و۱۳۸ ۱۱۳۹ ۱۱۸۳۱۷ 
۷ ۱۱/۰۸6۵ و۱۳۸ ۱۱۸۳۰۸۷ 
۸ ۱-۱( ۶-۰۱ ,+ ۱۳2۵ 
٩‏ ۱-۸-۰۹( ۷" ",+ 
۲۰ :۱۳,۳ ۱:۷ :۱۳,۳ 
میانگین و۱۳ ۱۱/۰۹۳۳۹۵ ۱۳۸ 
بهترین جواب ۱*۵ 1*۸۶ ۱-۱۹۳۲ 
بدترین جواب ۱۸۰۹ ۱9۵ ۱/۸۸۶۲ 
انحراف معیار ۱-۱ ۱( ۱7۹:۰۱۹۳,+" 
ضریب تغییرات 3 ی ۳ 


سال بیست و سوم شماره یک ۱۳۹۰ نشریه مهندسی عمران دانشگاه فردوسی مشهد 


۷۸ 


این مساله با استفاده از نرم افزار 1.160 (نسخه 
۸ مدل شده و جواب ۰/۷۲۱۰ برای آن حاصل شده 
است. همچنین. این مساله را جلالی [۱۷] با بکارگیری 
الگوریتم سیستم جامعه مورچه ها و معرفی الگوریتم 
«ي ۸0۹ مورد بررسی قرار داد و با بکارگیری ۱۵۰ 
مورچه و در ۵۰۰ تکرار مقدار ۰/۹4۵ را برای این مساله 
بدست آورد. در این مقاله. بهترین جواب بدست امده از 
بکارگیری الگوریتم ۷102۸5 بدون اعمال مکانیزم ۸ 
برابر با ۰/۹۳۱۰۲۵ می باشد. همچنین. بهترین. جواب 
پدست امده با استفاده از روش اول و دوم اعمال مک‌انیزم 
8 بترتیب برابر با ۰۷۵۰۸۶ و ۰/۸۱۱۸۳۲می باشد 
که نشان دهنده بهبود جوابه؛ بدلیل بکارگیری مکانیزم 
۶ می باشد. مقایسه نتایج بان دهنده: آبتت که با 
بکارگیری الگوریتم 2101۸5 بدون اعمال مکانیزم ۸ 
و با استفاده از روش اول و دوم اعمال مک‌انیزم 1۸1 
جوابهائی نزدیک به جواب بهینه مطلق حاصل شده است. 
این مساله یک مساله پیوسته می باشد اما الگوریتم جامعه 
مورچگان یک الگوریتم گسسته می باشد و همین موضوع 
علت تفاوت جواب‌ها با جواب حاصل از نرم افزار 
0 می باشد. همچنین» در حالت بدون استفاده از 
مکانیزم 1۸1 الگوریتم دچار همگرایی زود رس شده 
است؛ در حالیکه اعمال مکانیزم 19۸8 موجب ایجاد شک 
در مکانیزم جستجو و فرار از گیر افتادن جواب‌ها در نقطه 
بهینه محلی و حرکت به سمت جواب های مناسبتر 
(نزدیک به بهینه مطلق) می شود. 


مساله بهره برداری برفابی از مخحزن سد._ در این مساله 
تابع هدف به صورت حداقل سازی کمبود توان تولیدی 


نسبت به ظرفیت نصب یروگاه تعریف شده است. 


916 تج 


)۳۱( ۱ ۳ ۱ - ده 


- ۱۹ 


که در را بطه فوق. , ۴ توان تولیدی در دوره زمانی ] و 


نشریه مهندسی عمرال دانشگاه فردوسی مشهد 


بهره‌برداری پهینه از محازل سدها با استفاده از مکانیزم ی 


قظرفیت نصب یروگاه می باشد. 

علاوه بر قیود منظور شده برای مساله بهره برداری 
ساده. قیود دیگری نیز برای مساله بهره برداری برقابی از 
مخزن می بایست منظور گردد. به عنوان نمونه به منظور 
محاسبه توان تولیدی, علاوه بر میزان آب رها شده از 


مخزن,. میزان بار آب موثر بر توربین‌ها نیز می بایست 


تعریف گردد. بنابراین» مساله بهره برداری بهینه برقابی 5 
مخزن سد دارای قیود غیر خطی بوده و فضای جستجوی 
آن به صورت غیر محدب (000۷96) 08) می باشد. به 


نمود: 


٩ 21 (۳۲ 


1+ 


(۲۳) ۷۱۳۲۱۱۵۰۵۵۴۲۵0۲۶ 0۲) نج ع ظ 


ر( ,2-1 /( ,11 + 11) < .ظط 
(۲۵( (,۲)6 <, ۲۲ 
)۳۹( وک 
۳۷( تسیک کم سدع 


که روا بط فوق. ۲ وزن مخصوص آب؛ 1 بازده 
نیروگاه؛ ,0 بار آبی موثر بر نیروگاه در دوره زمانی ؛ 
۴ ضریب کارکرد نیروگاه؛ ,00 ضریب تبدیل حجم به 
دبی در دوره زمانی ؛ ,آتراز مخزن در دوره زمانی ] 
(تابعی از حجم مخزن در دوره زمانی ]) و ,1۷۷1 تراز 
پایاب نیروگاه می باشد. سایر پارامترها پیش از این 
تعریف شده است. 

اولین گام حل مساله: تعربف گراف مساله می 
باشد. متغیر تصمیم در نظرگرفته شده در این مساله مقدار 
آب رها شده از مخزن سد بوده و گراف این مساله نیز 
هماتيل ک اف تعریفت ی پا شاه یی دار ستاو 
از مخزن (شکل ۵) می باشد. 

بان با تایه ال تک ره تن 
برداری برقابی از مخزن سد دز مورد بررسی قرار گرفت. 
اطاتغات اوه ام سا نکن ات ان رکه هه اسست: 


سال بیست و سوم شماره یک. ۱۳۹۰ 
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همچنین. به منظور محاسبه میزان بار ابی موثر بر توربینها؛ 
با استفاده از اطلاعات حجم- ارتفاع مخزن. یک تابع حند 


جمله ای درجه سوم به صورت زیر بر آنها برازش داده 


شده است: 
٩‏ + 97 وج اجه عح 11 
,87205 0.05 < وه ,24983364 < 2 
۴ -< 4 ,1.37710 < 6 
)۳۸( 


برای اعمال قیود این مساله نیز از روش تابع جریمه 
ای که در مساله بهره برداری ساده شرح داده شد. استفاده 
شده است. مقدار ضریب جریمه در نظر گرفته شده در 
۳ 
منظور تعریف گراف مساله» محدوده مجاز مقدار آب رها 
شده از مخزن ([۱۰۰۰ ۰]) به ۱۱ دسته مساوی تقسیم 
شده است. مقدارمطلوب پارامترهای الگوریتم ۸5 
که از طریق آزمون و خطا بدست آمده اند. در جدول (۸) 
نشان داده شده اند. در اين مساله, به دلیل عدم امکان جدا 


سازی پارامترهای مقدار آب رها شده از مخزن سد و تراز 


۷۹ 


مخزن (که تابع حجم مخزن می باشد) در تابع هدف. 
امکان تعریف اطلاعات کاوشی برای هر یک از جزء 
مسیرها وجود نداشته و بنابراین ضریب تاثیر آن صفر 
منظور شده است (8-0). جدول )٩(‏ نشایج حاصل از 
حل این مساله بدون استفاده و با استفاده از روش دوم 
اعمال مکانیزم ۸۴ را به ازای ۱۰ اجرای مختلف, نشان 
می دهد؛ که در آن از طول بازه ۰/7 به منظور تعیین 
فضای جستجو جدید مساله در هر تکران استفاده شده 
1 
در نظر گرفته شده است. در این مساله از روش دوم 
اعمال مکانیزم ۶ سستفاده شده است. مقای‌سه نتایج 
نشان می دهد که اعمال مکانیزم 1۸ تاثیرقابل 
ملاحظه‌ای در بهبود جوابها داشته است. 


جدول ۸ مقادیر مطلوب پارامترهای الگوریتم ۸/۸5 


تعداد مورجه ها 


0 8 0 امه ۳ 
۱۰۰ ۱ ۰ ۰/۹۵ ۰/۹ 


جدول ۹ نتایج مساله‌بهره برداری برقابی از مخزن بدون استفاده و با استفاده از روش دوم اعمال مکانیزم ۱۱۳24 


تعداد اجرا 
۱ ۸۳/۹ 
۲ ۸۳۱ 
۳ ۸۳۲ 
3 وا 
۵ ۸/۳۷۷۹ 
1 ۸:۳۹ 
۷ 1-۹ 
۸ ۸-۸ 
۹ ۸-۸۳ 
۷۰ 2 
میانگین 1/۳۹۳ 
2 
۸۳۷۸۳ 
۳-۲( 


"۱۰۵ 


سال بیست و سوم شماره یک» ۱۳۹۰ 


بدون استفاده از مکانیزم ۱۱۳۹۹4۹ 


با استفاده از روش دوم مکانیزم ۱۳۹4 
۱-۹۶ 
۱/۵ 
کرک ۱/۸ 
۱۱۵ 
۱/۶ 
ات۱۸۵ 
۱/۳/۳۳ 
۱/۵۳۹ 
۱۹۲ 
۱/2 
۱2 
۱/۳/۳۳ 
1-۱ 
۱۹۹۹۹ 


"++-۳۳۵ 


تشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


نحوه همگرایی جواب‌ها برای بهترین جواب به 
دست آمده» در دو حالت بدون استفاده و با استفاده از 
روش دوم اعمال مکانیزم ۱۳4 در شکل (ع تیان داده 


شده است. همچنین. مقادیر انرژی تولیدی و آب رها 


استفاده از مک‌تیزم0۸6 0 
بدون استفاده از مکانیزم0۸۵ ---- 


1000 


تعداد تکرار 


بهره‌برداری بهینه از مخازن سدها يا استفاده از مکانیزم .. 


شده از مخزن در ماه‌های مختلف. در دو حالت بدون 
استفاده و با استفاده از روش دوم اعمال مکانیزم ۱۳۹4 به 
ترتیب در شکل‌های (۷) و (۸) نشان داده شده است. 


د 14 
13 
[ 12 
1 3 
3 
زاء 
1 دگِّ 
1 
[ و 
( 8 
7 
0 


شکل ٩‏ نحوه همگرایی جواب‌ها در مساله بهره برداری برقابی از مخزن سد دز در دو حالت بدون استفاده و با استفاده از روش دوم اعمال 


استفاده از مکنیزم0۸ --.6--- 
بدون استفاده از مکنیزم0(] هه 


40 50 60 


مکانیزم۸ 


30 
دوره زمانی(ماهانه) 


20 


انرژی تولیدی(۱۷۷) 


100 3 


0 10 


شکل ۷ مقادیر انرژی تولیدی در مساله بهره برداری برقابی از مخزن سد دز در دو حالت بدون استفاده و با استفاده از روش دوم اعمال مکانیزم 


نشریه مهندسی عمرال دانشگاه فردوسی مشعد 


9۳۹ 
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استفاده از مکنیزم 0۸6 -6.۰--- 
بدون استفاده از مکنیزم00۸6] هه 


40 50 60 


30 
دوره زمانی(ما هانه) 


۸ 


1000 


9 9 
فك فِ 
ت‌‌ِ ت‌‌ِ 


مقادیر رها شده از مخزن(۲61/0) 
جد 
8 


دح 
9 
ت‌‌ِ 


0 10 20 


شکل ۸ مقادیر آب رها شده از مخزن در مساله بهره برداری برقابی از مخزن سد دز در دو حالت بدون استفاده و با استفاده از روش دوم اعمال 


مکانیزم ۸ 


این مساله با استفاده از نرم افزار 1.1۱630 (نسخه 
۸ مدل شده و جواب ۷/۳۷۲ برای آن حاصل شده 
است. همچنین, این مساله را جلالی [۱۷] با بکارگیری 
الگوریتم سیستم جامعه مورچه ها و معرفی الگوریتم 
«ی۸05 مورد بررسی قرار داد و با بک‌ارگیری ۱۰۰ 
مورچه و در ۵۰۰ تکرار مقدار ۲۷/۲ را. که یک جواب 
ناشدنی است. برای این مساله به دست آورد. مقایسه 
نتایج جدول )٩(‏ نشان‌دهنده آنست که بهترین جواب به 
دست آمده از بکارگیری الگوریتم ۷1۳۸5 بدون استفاده 
از مکانیزم 7۸8 برابر با ۸/۱۰۶ می باشد. همچنین, 
بهترین جواب بدست امده از استفاده از روش دوم اعمال 
مکانیزم 2۸8 برابر با۷/۲۷۷ می باشد. که از جواب 
حاصل از مدل سازی در نرم افزار 113600 هم بهتر بوده 
است و نشان دهنده بهبود جواب‌ها به دلیل بکارگیری 
مکانیزم ۸1 می باشد. مقایسه نتایج بان دهفله انیس 
که با بکارگیری الگوریتم 1۷6۸5 بدون اعمال مکانیزم 
8 و استفاده از روش دوم اعمال مک‌انیزم ۸۴ 


جواب‌هائی نزدیک به جواب بهینه حاصل شتلاه آمنیتا: این 
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مساله یک مساله پیوسته می باشد اما الگوریتم جامعه 
مورچگان یک الگوریتم گسسته می باشد و همین موضوع 
مساله بهره برداری برقابی یک مساله غیر خحطی غیر 
محدب می باشد و لذا جواب بهینه حاصل شده از نرم 
امکان پیدا کردن جواب بهینه مطلق به دلیل مشکلات حل 
مسایل با استفاده از روشهای برنامه ریزی غیر حطیء غیر 
ممکن می باشد. همچنین در حالت بدون استفاده از 
مکانیزم ۸ الگوریتم دجار همگرایی زود رس شده 
است؛ در حالی‌که اعمال مکانیزم ۸ موجب ایجاد 
شک در مکانیزم جستجو و فرار از گیر افتادن جوابها در 
نقطه بهینه محلی و حرکت به سمت جواب های مناسبتر 
(نزدیک به بهینه) می شود. 


نتیجه گیری 


در این تحقیق مکانیزم مناسبی موسوم به تظریف تطبیقی 


تشریه مهندسی عمرال دانشگاه فردوسی مشعهد 


۸۲ 


جامعه مورچه ها برای حل مسائل بهینه سازی پیوسته 
ارائه شد. برای اعمال این مکانیزم دو روش پیشنهاد شد. 
نتایج این تحقیق نشان داد که بکارگیری مکانیزم تظریف 
تطبیقی قطعی در الگوریتم بهینه سازی جامعه مورچه ها 
در جهت نزدیکتر شدن به جواب های بهینه در حل 
مسائل بهینه سازی پیوسته. که یکی از چالشهای این 


الگوریتم محسوب می شود. بسیار موثر بوده انبتت ار 


بهره‌برداری پهینه از محازل سدها با استنفاده از مکانیزم ی 


قابلیت‌های مکانیزم پیشنهادی. با بکارگیری الگوریتم 
سیستم مورچه بيشینه - کمینه, در حل چند مساله ریاضی 
پیچیده و در مساله بهره برداری بهینه ساده و برقابی از 
مخزن سد دز مورد آزمون قرار گرفته و نتایم حاصل با 
نتایج حاصل از بکار گیری سایر روشها مقایسه شنک استت: 
مقایسه نتایج نشان‌دهنده بهیود نتایج در نتیجه اعمال 
مکانیزم پیشنهادی بوده است. 
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۷ جلالی محمد رضاء "طراحی و بهره برداری بهینه از هیدرو سیستمها با الگوریتم جامعه مورچه ها یک رهیافت فراکاوشی جدید" رساله 
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سال بیست و سوم شماره یک ۱۳۹۰ نشریه مهندسی عمران دانشگاه فردوسی مشهد 


۸ بهره‌برداری بهینه از مخازن سدها با استفاده از مکانیزم ... 
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